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STEADY-STATE COMPOSITION OF A LOW-DENSITY 
NONEQUILIBRIUM HYDROGEN PLASMA 
by Roman Krawec 
Lewis Research Center 

SUMMARY 

The steady- state composition of a hydrogen plasma is computed for both Maxwellian 
and non- Maxwellian electron energy distributions for the case in which only volume ef- 
fects need be considered. The electron energy distribution function is so defined that it 
can be smoothly varied from Maxwellian to Druyvesteyn. Steady-state plasma equations 

are derived and are solved numerically for selected forms of this energy distribution. 
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The results apply for plasmas that have a density of less than 10 particles per cubic 
centimeter. The gas is shown to be fully ionized for all distribution functions considered 
when the kinetic temperature of the electrons is above 4 electron volts. Certain regions 
of electron temperature are shown to exist in which Maxwellian and non-Maxwellian dis- 
tributions predict substantially different plasma compositions. 


INTRODUCTION 

The current interest in the production and properties of high-temperature plasmas 
has prompted renewed studies of the relation between temperature, energy distribution, 
and composition of such plasmas. One of the earlier attempts was made by Saha (refs. 

1 and 2), who formulated a relation between the temperature and the degree of ionization 
of a three-component (neutral atoms, ions, and electrons) plasma. Saha assumed the 
three components to behave as ideal monatomic gases and assumed complete thermo- 
dynamic equilibrium between the particles and the radiation field. An extension of Saha's 
work that allowed the radiation to escape from the plasma is described in reference 2. 

The analyses mentioned so far are applicable only to monatomic gases. A further 
treatment, with thermodynamic equilibrium assumed, was made by Gross and Eisen 
(ref. 3). This treatment considers a plasma composed of hydrogen atoms, molecules, 



atomic ions, negative ions, and electrons. The existence of molecular ions is not taken 
into account. 

Dewan (ref. 2) treated the case of a steady-state, nonequilibrium, three-component 
plasma and showed that the form of the electron energy distribution plays an important 
role in the relation between the electron temperature and the plasma composition. 

This report extends the work of Dewan to the case of a plasma composed of hydrogen 
molecules, atoms, molecular ions, atomic ions, and electrons. The method used is 
similar to that of Drawin (ref. 4). Both the effects of electron temperature and electron 
distribution are considered. The electron energy distribution function is defined in such 
a manner that it can be smoothly varied from Maxwellian to Druyvesteyn. 

The model used imposes certain restrictions on the applicability of the results. The 
first assumption is that the plasma composition is determined by processes within the 
plasma volume and that wall recombination is negligible. Since this assumption is not 
valid for many laboratory plasmas, the results of the present approach would not apply 
to those plasmas. 

A further restriction in the model used herein is that excited states are neglected. 
This restriction is valid if the probability of a collision of an electron with an excited 
molecule or atom is small compared with the probability of a collision with a ground- 
state molecule or atom, which is the case if the plasma density is not too high or if the 
plasma is optically thin. 

The calculations are performed for molecular hydrogen densities of 10® and 10^^ 
particles per cubic centimeter and kinetic temperatures from 0. 4 to 6. 0 electron volts. 
Experimental values have been used for the cross sections and collision frequencies 
whenever possible. 


SYMBOLS 


A 

a 

B 

b 

C, Cj, C£ 
D 
E 
F 


rate coefficient, defined in appendix 
normalizing factor, defined in eq. (15) 
rate coefficient, defined in appendix 
defined in eq. (15) 

rate coefficients, defined in appendix 
rate coefficient, defined in appendix 
rate coefficient, defined in appendix 
rate coefficient, defined in appendix 
reaction rate of process (fl) 
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Fg reaction rate for process (f2) 

k Boltzmann constant, ergs/°K 

m electron mass, g 

3 

n number of particles per cm 

3 

n m original number of particles per cm 

T temperature, °K 

t time, sec 

v electron velocity, cm/sec 

a recombination coefficient 

(a!j> radiative recombination coefficient 

(a g) collisional recombination coefficient 

j3 3. 5xl0 16 /(torr)(cm 3 ) 

e electron energy, ergs 

(e) average energy, ergs 

a cross section, cm 3 

r kinetic temperature, defined in eq. (16) 

v collision frequency, (torr”*)(sec~^) 

* vn mo^ 

Subscripts: 

ac atomic collisional 

ai atomic ion 

ao neutral atom 

ar atomic radiative 

dr dissociative recombination 

e electron 

id molecular ion dissociation 

mi molecular ion 

mo neutral molecule 

rc radiative collisional 
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THEORY 

Model 


As previously mentioned, the plasma is considered to be optically thin and suf- 
ficiently large or well contained so that wall effects are negligible. 

Only the following reactions are considered: 

(a) H + e — H + + 2e (Atomic ionization) 

(b) H 2 + e — H 2 + + 2e (Molecular ionization) 

(cl) e + H + — H (Radiative recombination) 

(c2) e + e + H + — H + e (Collisional recombination) 

(d) e + H 2 + — H + H (Dissociative recombination) 

(e) H 2 + + e— H + +H + e (Dissociation of molecular ions) 

(f 1) H + H + Hp — H« + Hq') 

V (Three- body atomic recombination) 

(f 2) H + H + H — H 2 + H J 

A search of the literature indicated that all other reactions were negligible com- 
pared with those included with the possible exception of ionization of excited atoms or 
molecules. The neglect of this type of ionizing collision can be justified if the density to 
which the model is applied is properly restricted. Since the decay time of an excited 
state is of the order of 10“ 8 second (ref. 5, p. 9), it is sufficient to impose the condition 
that the time between ionizing collisions be large compared with 10“ 8 second. With the 

use of the maximum value of the ionization frequency for the hydrogen molecule, the 
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model is suitable for densities less than about 10 particles per cubic centimeter. For 

_6 

these densities, the time between ionizing collisions is at least 10 second. 

At the lower temperature and higher density end of the region treated herein, the 
plasma may tend to become optically thick for the resonance lines. Thus, the assump- 
tion of optical thinness must be carefully checked by independent methods. Other as- 
sumptions are that the plasma is considered to be in a steady state, that the ions and the 
neutrals are at room temperature (charge exchange reactions can be neglected), and that 
the original gas is composed of diatomic molecules. 


PRODUCTION RATES 

The number of atomic ions formed per unit time per unit volume by process (a) is 


^< n al> = n ao / “al v dn e * 


an n A 
e ao 


(i) 
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where a is a normalizing constant that depends on the distribution function, and A is 
the rate coefficient for process (a). (All rate coefficients are evaluated in the appendix. ) 

Similarly, the number of molecular ions formed per unit time per unit volume by 
process (b) is 


J t ( Vi> = 




( 2 ) 


The loss rate of atomic ions due to radiative recombination (process (cl)) is 


--(n a r) = <«i> n ai n e san e n ai C l 


(3a) 


An integration is not performed for this case since only the total recombination co- 
efficient was known. The cross section for this process, however, is expected to ap- 
proach infinity as the energy approaches zero and is a monotonically decreasing function 
of energy. Process (d) also exhibits this behavior, and the appendix shows that the form 
of the distribution function has little effect on the result. Thus, using a rate coefficient 
that is independent of the energy distribution is not believed to introduce a large error. 

A similar consideration holds for collisional recombination (process (C2)), which is 
written 


^ n ac> = <“2> n ai n e s an ai n e C 2 < 3b ) 

The rate of atom production due to dissociative recombination (process (d)) is 

J t (n dr) = 2n mi r CTdrV ^ " 2an e n mi D < 4 > 

J Q 

The rate of production of atoms (or atomic ions) by dissociation of molecular ions 
(process (e)) is 


d 

— (n. = n / a.jV dn s an n -E 

v id' mi I id e e mi 

JfX 


( 5 ) 


The rate of loss of atoms due to processes (f 1) and (f2) is 
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( 6 ) 


^ ( n ao^ ^ n ao I1 mo F l + 2n ao F 2 2n ao F 


Here, and F 2 are the reaction rates for processes (fl) and (f2), respectively. 


Genera! Steady-State Equations 

At steady-state conditions, the rate of atomic ion production (processes (a) and (e)) 
is equal to the rate of atomic ion consumption (processes (c) and (c2)). Thus, 

n ao^ + n mi E “ n ai^ ^ ‘ 


where 


C - Cj + C 2 n e 


The reactions that give rise to an increase of electrons (processes (a) and (b)) must 
be balanced by the reactions in which an electron is captured (processes (c) and (d)). 
Thus, 

n A+n B = n .C + n .D (8) 

ao mo ai mi v } 

The conservation of neutral atoms requires that 

2n 2 F 

on 

n ao A+ — = n ai C+ “ml (2DtE) < 9 > 

an e 

The only other equations that are needed are those specifying conservation of charge 
and conservation of mass: 


n e = n ai + Vi 


V * Vo + Vi + j (n ao + n ai> 


Combining equations (7) and (8) results in 


irbuvtt!” on 


mi 

mo 


B 

D + E 


Combining equations (7) and (9) gives 


n 2 F 
ao 

an 


= n mi( D + E > 


( 12 ) 


(13) 


Equations (7), (10), (11), (12), and (13) are the general steady-state equations that 
must be solved simultaneously for the amounts of the various components. 


Energy Distribution Function 

In order to evaluate the coefficients needed to solve the general steady- state equa- 
tions, the following expression for the electron energy distribution was arbitrarily 
chosen: 


dn oane-^V^de 


(14) 


where 1 < b < 2, and a is a normalizing factor given by 


a = 


/ e-<*/ T >V/2 d£ r 3/2 rg\ 


(15) 


The factor b allows the energy distribution to be varied from Maxwellian (b = 1) to 
Druyvesteyn (b = 2). The electron kinetic temperature r is related to the electron tem- 
perature T in the following manner: 


3 kT = ( e ) = a f e 3/2 e -(e/T) b de 

O •'A 


(16) 


Combining equations (15) and (16) produces 
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Normalized energy distribution function, erg" 1 Ratio of Maxwellian temperature to kinetic 

e 08 temperature, kT/x 


kT = 


(17) 




Figure 1. - Effect of energy parameter b on ratio 
of Maxwellian temperature to kinetic tempera- 
ture. 



Figure 2. - Energy distribution function. Temperature kT, 
1 electron volt. 


This method defines an electron kinetic 
temperature that is related to the average 
energy of the electrons in the same manner 
as for a Maxwellian distribution. This rela- 
tion is given in figure 1 for the range of 
values of b considered herein, while fig- 
ure 2 shows the general shape of the distri- 
bution function plotted for two values of b 
at a value of kT of 1 electron volt. 

RESULTS AND DISCUSSION 

The rate coefficients were evaluated for 
different forms of the electron energy distri- 
bution function by numerical integration on an 
IBM 7094 computer. The values obtained for 
these rate coefficients were then used to 
solve the general steady- state equations for 
the fractional amounts of the various parti- 
cles present. 

The fractional amount of the particles 
present is presented in figures 3 and 4 for 
different energy distributions (different 
values of b) as a function of electron tem- 
perature. The composition for a Maxwellian 
energy distribution is given for initial densi- 

Q IQ 

ties of 10 and 10 particles per cubic centi- 
meter in figure 3. The higher density case 
(fig. 3(b)) was solved both with and without 
the inclusion of collisional recombination. 

This term has little effect in the final results. 
Neglecting this process had the advantage of 
simplif ying the computer program consider- 
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(a) Original number of particles per cubic centimeter, 10 9 . (b) Original number of particles per cubic centimeter, 10 13 . 

Figure 3. - Composition of hydrogen plasma for Maxwellian energy distribution. Energy parameter, b - 1 

ably and was done for all the lower density cases (figs. 3(a) and 4). 

The effect of changing the form of the electron energy distribution is illustrated in 
figure 4, in which the fractional amounts of the different particles are given as a function 
of electron temperatures for four energy distributions. Figure 4(c) shows that a 
Maxwellian plasma would be fully ionized for a temperature of 2. 5 electron volts, while 
the fractional amount of atomic ions for a Druyvesteyn plasma would be 0. 57. The latter 
plasma would be completely ionized at temperatures above 3. 85 electron volts. 

Figure 4 also shows that the type of energy distribution function does not change the 
shape of a given curve to any extent but merely causes a shift of the curve along the tem- 
perature axis. For a given temperature, this figure also shows that the fractional 
amounts of a given particle can change by orders of magnitude as the distribution function 
varies. The region where most changes in ionization take place lies between 0. 5 and 3. 0 
electron volts. In this region, the energy distribution of the electrons is as important as 
the electron temperature in determining the degree of ionization. 

In practice, it is doubtful whether the energy distribution will be sufficiently well 
known to be of use in determining percentage of ionization. With an independent mea- 
surement of the percentage of ionization and the electron temperature, however, an 
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Temperature, kT, eV 

<a> Molecular neutrals. (b) Atomic neutrals. ( C ) Atomic ions. 

Figure 4. - Fractional amount of molecular neutrals, atomic neutrals, and atomic ions for different forms of distribution function. Original number of particles per cubic 
centimeter, 10 y . 


TABLE I. - EFFECT OF CHANGES IN RATE COEFFICIENTS ON REPORT RESULTS 


[Original number of particles per cubic centimeter n m , 10 ; energy parameter b, 1; rate coefficient Cg, 0. ] 


Kinetic 

temperature, 

t, 

eV 


Rate coefficients 


Fractional amounts 


A, 


B, 


Ci, 


D, 


E, 


-1 


-1 


-1 


-1 


C 

cm /sec 


cm® /sec 


Molecular Atomic Atomic 


neutrals, 

n mc/ n m 

neutrals, 

n /n 
ao' m 

ions 

n ai/ n m 

0.485 

1. 03 

9. 60X10' 5 

.782 

.415 

3. 87x10'® 

. 137 

1.73 

1. 62x10' 3 

.485 

1. 03 

9. 59x10' 5 


0.8 


7. 940x10 


■ 35 


2. 780X10 


-34 


8. 512X10 
8. 512x10 
8. 512x10 
8. 512x10 


-31 

30 
32 

31 


2. 84x10 
2. 84x10 
2. 84x10 
2. 84x10 
2. 84x10 
2. 84X10 


7. 005x10 


•32 


3. 542x10” 50 1.188X10 


-47 


7. 005X10 
7. 005x10 
7. 005x10 


31 
-33 

32 


3. 542X10 
3. 542x10' 
3. 542x10 
3. 542x10' 


-49 

,-51 

-50 

50 


1. 188X10 
1. 188X10 


-46 

48 


.487 

.485 

.789 

.137 


.415 

1.73 


1.46X10' 3 

1.99 

1.43x10" 2 

1. 97 

1. 38x10' 4 

1.93 

1.46X10' 3 

1 

1.99 

i 


9. 54x10 
9. 56X10 


3. 88X10 
1. 61x10 


-5 

-5 

-5 

-4 


1.0 


3. 72x10 


-33 


1. 286X10 


-32 


1. 064x10 
1. 064x10 
1. 064x10 
1. 064x10 


-30 

29 
-31 

30 


3. 55x10 
3. 55x10' 
3. 55x10 
3. 55x10 
3. 55X10 
3. 55X10 


-20 

,-20 

,-20 

-19 

-21 

20 


1. 976X10' 


-30 


4. 95x10' 


-50 


1. 66X10 


■47 


1. 976x10 
1. 976x10 
1. 976X10 


-29 
•31 
■ 30 


4. 95x10 
4.95x10 
4.95x10 
4. 95x10' 


,-49 

-51 

-50 

-50 


1. 66x10 
1. 66x10 


•46 

-48 


1. 42x10 
1. 47x10 


1.97 

1.99 


6. 94x10 
6. 90x10 
6. 74x10 
6. 94x10 
6. 92x10 
6. 94x10 
6. 92x10 
6. 92X10 
6. 92x10 
6. 86x10 
6. 95x10 


1.5 


5. 77x10" 31 2.864x10 


30 


1. 596x10 
1. 596x10 
1. 596x10 
1. 596X10 


5. 325x10 
5. 325x10 
5. 325X10 
5. 325X10 
5. 325X10 
5. 325x10 


20 

,-20 

-20 

-19 

21 

26 


1. 918x10 


t 

1.918x10’ 
1.918x10’ 
1. 918x10“ 


-28 


9. 094x10 


-50 


3. 05x10 


-47 


9. 094x10 
9. 094x10 
9. 094x10 
9. 094x10 


,-49 

-5! 

,-50 

-50 


6.36X10" 8 1.47 
1. lOxlO' 6 1.93 

5. 52x10” 10 4. 32X10" 1 

6. 36X10" 8 1.47 


3. 05x10 
3. 05x10 


6. 36x10” 
6. 36X10 


-9 


5. 32X10 
7. 00x10 
1. 57 
5. 30x10 
5. 32x10 
5. 32x10 
5. 30X10 
5. 32X10 



’ 'effective” b can be obtained to 

describe the electron distribution. 

A comparison of the present work 

(fig. 5) with the Saha equation 
13 

(n = 10 ) and with the results of 

Gross and Eisen (n = 2. 7x10^; 
ref. 3) shows that a higher tempera- 
ture is needed to obtain complete 
ionization than is predicted from their 
thermodynamic approach. An exact 
comparison cannot be made with 
reference 4 since the density was not 
specified. 

The model used herein is be- 
lieved to be a valid approximation for 
a well-confined plasma that is formed 
by nonthermal means. Such a plasma 
could be formed by a radio frequency 
electrodeless discharge in a strong 
magnetic field. For such a case, the 
use of an equation derived for a 
plasma in thermodynamic equilibrium 
could easily lead to extremely large errors in predicted plasma composition. 

An idea of the sensitivity of the method to the values of the various cross sections 
used was obtained by changing the values of the rate coefficients used by factors of 0. 1 
and 10. These data are presented in table I for different values of the kinetic tempera- 
ture. The first row for each of the three different values of t gives the results that 
were presented in figure 3(a). The other rows show the effect on the results when various 
rate coefficients are changed. Radiative (cl) and three-body (f2) recombinations are seen 
to be the most important. 



Figure 5. - Comparison of fractional amounts of atomic particles from 
three different calculations. Number of particles per cubic centi- 
meter, n m , 10 . 


CONCLUDING REMARKS 

The composition of a plasma in which wall processes can be neglected was calcu- 
lated for various nonthermal energy distributions. The method consisted of deriving a 
set of steady- state equations that depended on the electron energy distribution function 
along with appropriate cross sections or collision frequencies. 
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The reactions that were found to be important in determining the composition of the 
plasma were 

1. Electron bombardment of atoms to form atomic ions 

2. Electron bombardment of molecules to form molecular ions 

3. Radiative recombination of atomic ions 

4. Dissociative recombination of molecular ions 

5. Dissociation of molecular ions 

6. Three-body recombination of atoms to form molecules 

The temperature needed to dissociate and ionize the gas completely increased as the 
distribution was made more non- Maxwellian. A prediction of plasma composition for a 
given temperature could easily differ by large values if the plasma were non- Maxwellian 
instead of Maxwellian. The type of energy distribution did not change the shape of any of 
the curves but only shifted them toward higher temperatures. 

All the distribution functions considered predicted complete ionization when the tem- 
perature was above 3. 85 electron volts. 

An equation derived from a thermodynamic approach can lead to large errors in pre- 
dicting the composition of a nonthermal plasma. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 21, 1966. 



APPENDIX - EVALUATION OF RATE COEFFICIENTS 


From equation (1), 


A 



CaiVdn 


e 


(Al) 


The value of a. as a function of energy was taken from reference 6 (p. 434) and is 

SIX 

given in figure 6. Using this cross section along with equations (14) and (15) for dn g 
and a allowed the coefficient to be evaluated. 

From equation (2) 


B s 


1 

a ^ n e 



(A2) 


Values of t' mo were obtained from reference 7 and are plotted in figure 7. 

The recombination of an electron with an atomic ion can proceed either by radiative 
or collisional recombination processes (cl) and (c2). These two cases are treated sepa- 
rately. 

Little experimental work has been done on atomic radiative recombination coeffi- 
cients. The experimental data of references 8 and 9 appear to justify the use of the 
theoretical formulation given by Dewan (ref. 2), wherein 

1. 6xl0" 32 

"" = <e> 

Use of this equation gives 
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<«!> 

a v or 
ar ar 

f2(e>" 

ly/2 7. 51x10" 19 r „/ 3 \ 

a 

a a 

L mJ 

b r u 


(A3) 


A combined radiative -collisional recombination coefficient is treated in refer- 
ence 10. * This combined coefficient is shown to approach a pure radiative recombina- 
tion coefficient at low densities and become collision dominated at higher densities. The 
data presented were used to infer that Q! rc » + <*2 n e 


6 

E 

!> 


o 

’«/» 




<« 2 > = 


1. 27X10 

<€> 1 * * * 5 * 


•86 


Hence, 



1. 27x10" 86 r (^j 

br 7 / 2 


(A4) 


From equation (4) 




a . v dn. 
dr e 


(A5) 


The theoretical calculation of the re- 
combination coefficient for dissociative re- 
combination was made by Bates (ref. 11) 
who found that 


Figure 7. - Collision frequency for ionization of hydrogen 
by electron bombardment (ref. 7). 


"dr = 


C onstant 

. 1/2 


(A6) 


1 A three -body recombination coefficient could have been derived independently of 

Bates (ref. 10) by using the principle of detailed balancing. (See, e. g. , J. B. Hasted, 

Physics of Atom Collisions, Butterworth & Co. , 1964, or Hans R. Griem, Plasma 
Spectroscopy, McGraw-Hill, Pub. Co., 1964.) However, the method used was deemed 

sufficient for the purposes of this report. It may be objected that the use of a recombina- 

tion coefficient derived from Bates is inconsistent with the neglect of excited states in 

setting up the rate equations. Strickly speaking, this is true. Since the expressions are 
derived for an optically thin plasma, however, the effect of excited states should be neg- 

ligible. 
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H 

An order -of -magnitude -type calculation of the constant gave Bates a value of a dr = 10 
cubic centimeter per ion second for an electron temperature 250° K. Biodi and Brown 

a 

(ref. 12) found an experimental value of a^ r = 2. 5x10 cubic centimeter per ion second 
for the electron -ion recombination coefficient at a temperature of 300° K but were not 
sure of their electron loss mechanism. The cross-sectional data of Rose and Clark 
(ref. 13) gave a value of the recombination coefficient near that given in reference 11. 
Their value along with equation (A6) gives 

< € / r)b de = 2. 22x10" 8 ^ T^j 

This rate coefficient is not strongly dependent on the form of the energy distribution func- 
tion. From equation (5), 


/ OO 

« 


, dn = Constant x 
dr e 


f 



an 


e 



a id v dn 


e 


(A7) 


Experimental values of <7^ could not be found in the literature. Ivash (ref. 14) per- 
formed a theoretical calculation of oj d and included the effects of exchange interactions. 
These interactions affect the cross sections only at energies below 50 electron volts, 
while for higher energies, cr^ may be expressed by 


"id" 


10.75 


0.917 

i- 12 


(A8) 


1. 6X10 

Ivash calculated separate cross sections when a prior interaction or postinteraction 
was assumed. The values found by assuming postinteraction seemed too high from energy 

considerations; hence, the prior interaction 
values were used for energies less than 71 
electron volts (fig. 8), while equation (A8) 
was used for higher energies. 

From equation (6), 



F = F.n^ + F 0 n „ 
1 mo 2 ao 


(A9) 


80 


Figure 8. - Cross section for dissociation of molecular ions by 
electron bombardment (ref. Id). 


Reference 15 gives reaction rates for 
the two three -body processes considered. 
The values from reference 15 give 


F = 2. 76X10 


32 n mn + 9 - 25xl0 ’ 3 ° n n n 

mo ao 
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